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The  distortion  due  to  the  turbulent  boundary  layer  becomes  magnified  by  reduction 
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AERO-OPTICS  AT  SHORTER  WAVELENGTHS 

by 

Distinguished  Professor  Allen  E.  Fuhs 

Department  of  Aeronautics 

Department  of  Physics  and  Chemistry 

Naval  Postgraduate  School 

Monterey,  California,  93940 

I .  INTRODUCTION 

An  extensive  technological  base  for  aero-optics  has  been  developed  by  the 
laser  community.   The  research  and  development  work  was  oriented  to  the  C0~ 
laser  with  wavelength  of  10.6  microns.   In  recent  times,  lasers  operating  at 
much  shorter  wavelengths,  e.g.,  iodine  laser  at  1.315  microns  [1]  or  the 
excimer  lasers  in  the  blue  end  of  the  visible  spectrum  [2],  have  shown  potential 
for  superior  performance.   The  information  concerning  laser  turrets  on  aircraft 
needs  to  be  reexamined  for  shorter  wavelengths. 

Section  II  discusses  the  far-field  intensity  using  an  equation  which  was 
developed  by  Holmes  and  Avizonis  [3].   The  equation  highlights  the  impact  of  a 
change  in  wavelength.   Section  III  discusses  scaling  of  turret  weight,  volume,  and 
drag.   Section  IV  discusses  the  scaling  of  phase  distortion  due  to  inviscid 
flow.   Section  V  considers  phase  distortion  to  aircraft  generated  turbulence. 
Finally  Section  VI  summarizes  the  paper  and  states  conclusions  relative  to 
unknowns  or  areas  of  uncertainty. 

II.  FAR-FIELD  INTENSITY 

An  equation  for  far-field  intensity  has  been  developed  by  Holmes  and 
Avizonis  [3].   The  equation,  which  is  useful  for  systems  analysis,  is 


2    RmT  T,T  K  exp[-(27T6    /A)2] 
.ttD  P  .   w  d  a    v  rms ,,v 

ff  '  '  *  2   2)  2  2 

4*  Z      (_L.)   [i  +  (5£D}  , 
lFDQ;   L±    K3X    }    J 

One  additional  term,  D/FD  ,  has  been  added  by  Fuhs  [4]  to  account  for  turbulence. 
The  meaning  of  each  term  will  be  discussed  now. 

The  first  term  is  the  far-field  intensity  on  axis  for  an  Airy  diffraction 
pattern.   The  mirror  diameter  is  D;  the  laser  power  at  the  output  aperture, 
i.e.,  laser  telescope  primary  mirror,  is  P;  the  wavelength  of  the  laser  is  A; 
and  the  range  to  the  focal  plane  is  Z.   The  shorter  wavelength  increases  lff; 
however,  one  may  want  to  decrease  mirror  size,  thereby  decreasing  turret  weight 
and  volume.   Further,  a  decrease  in  mirror  diameter,  D,  will  decrease  aerodynamic 
drag.   Consequently,  the  ratio  of  D  to  A  will  be  assumed  to  be  constant,  i.e., 

S  =  A/D  =  constant  (2) 

S   is  a  scaling  factor.   The  turret  will  have  a  characteristic  dimension  L.   An 


additional  assumption  is 


S   =  L/D  =  constant  (3) 


where  S   is  another  scaling  parameter. 
The  term  for  jitter  is 

1  +  (ff>2  =  1  +  <%>2  <« 

When  S   is  constant,  the  jitter  term  depends  on  the  line-of -sight  rms  angular 
beam  jitter,  a.   Many  different  phenomena  influence  the  magnitude  of  a  including 
the  following: 

aerodynamic  inputs  to  airframe 

aerodynamic  input  directly  to  the  laser  turret 

aerodynamic  input  due  to  control  surface  deflections 

unsteady  aerodynamics  in  vicinity  of  laser  turret 


unsteady  aerodynamics  due  to  open  port 

vibration  due  to  aircraft  components 

electrical  noise  within  servo  system 

laser  turret  mechanical  jitter 

interaction  of  turret  aerodynamics  with  remainder  of  aircraft 
Each  item  will  not  be  discussed.   When  the  laser  turret  shrinks  in  size  to  main- 
tain S1  constant,  the  value  of  a  must  be  maintained  constant  to  avoid  increased 
detrimental  effects. 

From  the  causes  of  jitter  listed  above,  the  problem  due  to  open  port  may  be 
eliminated  or  decreased.   With  another  wavelength,  i.e.,  1.3  microns  instead  of 
10.6  microns,  and  a  much  smaller  size,  a  window  may  become  feasible. 

The  quantity  R  is  the  loss  due  to  the  reflectivity  of  m-mirrors  in  the 

optical  train.   T  is  the  transmission  of  the  window  mentioned  in  the  preceding 
w 

paragraph.   The  quantity  T   is  diffractive  spillage  around  mirrors  in  the 

optical  train.   T,  is  a  function  of  Fresnel  number 
a 

2    2 
r    D      D 

m    __E  _  m  _   in  ,_. 

F   Xi       4AJI   4S  £  K   } 

where  r  is  the  radius  of  mirrors  in  the  optical  train,  and  I   is  the  spacing 
m 

between  two  mirrors.   If  one  scales  the  primary  mirror  by  ratio  S  ,  then  the 

diameter  of  a  mirror  in  the  optical  train,  D  ,  undoubtedly  will  be  scaled  to 

m 

a  smaller  size.   The  consequence  is  increased  spillage,  i.e.,  a  smaller  value 

for  T,.   One  assumes  that  I   remains  fixed, 
d 

The  transmission  of  the  atmosphere  is  T  and  is  given  by 

T  =  e"  YZ  (6) 


where  y  is  an  extinction  coefficient.   The  extinction  coefficient  is  the  sum 

of   four   terms 

Y=a+a+k+k  (7) 

a  s  a  s 

where  a      is  molecular  absorption,  a   is  molecular  scattering,  k   is  aerosol 
a  s  a 

absorption,  and  k   is  aerosol  scattering.   Each  of  the  terms  in  equation  (7)  is 

dependent  on  wavelength. 

2 
The  factor  K  is  a  consequence  of  nonuniform  intensity,  I,  watts/cm  ,  at 

the  primary  mirror  of  the  laser  turret.   If  S   is  fixed,  the  intensity  in  the 

far-field  does  not  change;  see  Appendix  A. 

Aero-optics  plays  an  important  role  in  the  factor  6   /A.   Of  course,  6 

r     v     J  v  rms  rms 

is  the  result  of  all  optical  aberrations  starting  within  the  cavity  of  the  laser 

device.   The  aerodynamic  window,  the  mirrors  in  the  optical  train,  and  the 

primary  telescope  mirror  affect  6    .   The  beam  propagates  through  a  turbulent 

boundary  layer  and  the  variable  density  field  generated  by  the  inviscid  flow  over 

the  turret.   Both  external  flows  contribute  to  6 

rms 

Obviously,  when  A  is  decreased,  6    must  be  decreased  by  a  like  amount. 

rms 

Otherwise,  Iff  is  reduced.   When  A  is  reduced  from  10.6  microns  to  1.3  microns, 

i.e.,  by  a  factor  of  8,  6    must  be  reduced  also  by  a  factor  of  8.   The  most 

rms 


difficult  aspect  of  scaling  to  shorter  wavelengths  will  be  to  reduce  6 


rms 


A  connection  exists  between  ip CC »H )  in  equation  (A-l)  and  6   /A  in  equation  (1)  . 

For  more  details  concerning  the  relationship  between  \b   and  6    ,  see  Born  and 

rms 

Wolfe  [5]. 

A  term,  D/FD  ,  appears  in  equation  (1)  to  account  for  turbulence.   The  term  is 
adequate  for  systems  analysis  but  is  not  appropriate  for  a  detailed  examination  of 
aero-optics.   The  question  of  turbulent  flow  is  deferred  to  Section  IV. 


III.   SCALING  TURRET  WEIGHT,  VOLUME,  AND  DRAG 

Turret  weight,  W,  can  be  scaled  using  the  equation 

w^  "  ^  (8> 

where  n  is  a  number  depending  on  the  turret  optics  and  D  is  output  aperture. 
Table  I  gives  the  value  of  n,  W  ,  and  D..  for  three  turrets. 


Table  I.   Parameters  for  Turret  Weight 


Type  of  Turret 


n 

W  ,  pounds 

D  ,  cm 

2.50 

1637 

60 

1.93 

1788 

60 

2.40 

2161 

60 

on-gimbal  Cassegrain 

heliostat 

coelostat 


Changing  from  a  CO  laser  to  an  iodine  laser  gives  a  wavelength  ratio  of  1/8. 
Assuming  S   is  fixed,  the  value  of  D  can  be  reduced  by  a  factor  of  1/8;  mirror 
diameter  can  be  reduced  from  60  cm  to  7.5  cm.   Use  the  numbers  in  equation  (8)  for 
an  on-gimbal  Cassegrain  telescope.   The  result  is 

1  2'5 
W2  =  1637  [-]    =9.04 

The  range  of  validity  of  the  scaling  law  of  equation  (8)  may  have  been  exceeded; 
even  so,  the  qualitative  conclusion  is  that  turret  weight  is  very  significantly 
reduced. 

Turret  volume,  V,  follows  a  scaling  law  of 

3 

V     D 

—  =  (— )  (9) 

Once  again,  significant  reductions  occur. 


Aerodynamic  forces  are  given  by  an  equation  of  the  form 


F  =  ^pu2ACF  (10) 

where  C  is  a  force  coefficient.   The  drag  is 

D  =  fpu2ACD  (ID 

where  C   is  the  drag  coefficient.   In  equation  (11),  A  is  the  cross  sectional 
area  of  the  turret  exposed  to  the  flow.   Hence  the  drag  varies  as  mirror  diameter 
squared. 

When  laser  turrets  weigh  in  tons,  only  one  turret  can  be  used  on  an  aircraft 
However,  when  the  turret  is  reduced  in  weight  and  volume  by  a  factor  of  100 
or  so,  many  turrets  can  be  used.   Instead  of  attempting  to  cover  27T-steradians 
with  one  turret,  each  of  the  small  turrets  can  cover  a  much  smaller  solid  angle. 
Furthermore,  the  turrets  can  be  located  at  many  different  places  on  the  aircraft. 
In  the  enthusiasm  to  add  turrets,  consideration  must  be  given  to  fire  control, 
pointing,  and  tracking. 

IV.   SCALING  OPTICAL  DISTORTION  DUE  TO  INVISCID  FLOW 

Three  lengths  are  of  relevance  to  the  scaling  of  optical  distortion  due  to 
inviscid  flow;  these  are  beam  diameter,  D;  wavelength,  X;  and  characteristic 
turret  dimension,  L.   One  length  is  selected  as  a  reference  length  to  form 
ratios.   Equations  (2)  and  (3)  indicate  D  is  used  as  the  reference. 

Compressible  flow  over  a  laser  turret  generates  a  spatially  varying  density 
field.   A  relation  exists  between  refractive  index  and  density 

n  =  1  +  (Kp/pQ)  (12) 


The   fractional   changes  are 


An  Ap 

n~=Kp  d3) 


By  solving  the  equations  of  motion  for  flow  over  a  turret,  the  density  field  can 
be  obtained.   Knowing  the  density  field,  the  refractive  index  can  be  calculated. 
The  optical  path  length,  OPL,  is  defined  as 


OPL  = 


=  /n(x,y,z)ds  (14) 


where  s  is  distance  along  a  ray.   The  optical  path  length  is  a  geometrical 
length  weighted  by  the  local  refractive  index.   The  change  in  OPL,  known  as 
optical  path  difference,  OPD,  is 

6  =  OPD  =  /n.(x,y,z)ds.  -  /n  (x,y,z)ds  (15) 

where  subscripts  i  and  r  indicate  the  i-th  ray  and  a  reference  ray.   The  quantity  6 
is  equal  to  the  optical  path  difference. 
Define  nondimensional  variables 

x  =  x'D 

y  =  y'D 

z  =  z'D 

s  =  s'D 
Substitution  of  the  nondimensional  variables  into  equation  (14)  gives 

OPL  =  D  f nU'.y'.z^ds'  (16) 

Likewise  the  optical  path  difference  becomes 

6  =  D  [n.ds!  -  D  fnds*  (17) 

J      i   i     J     r   r 

Hence  OPL  and  6  scale  as  mirror  diameter  D.   One  must  show  that  the  integrals 
in  equations  (16)  and  (17)  are  constant  when  turret  size  varies  so  that  S,  and  S 
are  constant. 


For  blunt  turrets,  the  equation  for  the  velocity  potential  in  compressible 

flow  is 

2  2 

1     9$     9(j)  9   9     _     9  g>  (18) 

a     9x^    9x,    9xJ3x,         9x.9x. 
i       Tc        x     k  j      j 

Introduce  nondimensional  variables 

9  =  9'u  D 


a  =  a'a 

00 


into   equation   (14).      The  result   is 


jgJL  3£1  9^1     924>'  924>' 

,2  " 

OA  .       OA.         OA  .  OA,  O/ 

j  j 


a'      9x!    9x/    9x|9x/ 


(19) 


A  length  scale  parameter  does  not  appear  in  equation  (19) .   Although  not  demonstrated 
here,  the  boundary  conditions  can  be  nondimensionalized  without  a  residual  length 
scale.   Consequently,  the  solution  <() '  is  independent  of  turret  scale. 
The  local  Mach  number  is  related  to  freestream  Mach  number  by 

u'u 

u        °° 

M  =  -  =  — : —  =  M  M'  (20) 

a   a  a     °° 

oo 

The  local  Mach  number  determines  the  local  refractive  index.   Since  M'  is 
independent  of  turret  scale,  the  identical  value  of  refractive  index  exists  at 
corresponding  points  in  two  flows  with  different  scale.   Consequently,  the  integrals 
in  equations  (16)  and  (17)  are  identical  for  the  two  flows  at  different  scale. 
The  phase  distortion,  also  known  as  optical  path  difference,  relative  to 
wavelength  is 


=  f/nidSl  -f/nrdSr  (21) 


Since  the  integrals  are  independent  of  scale,  and  since  S   is  fixed,  the  phase 

distortion  does  not  vary.   Likewise  6   /A  does  not  vary  with  turret  scale 

rms  J 

provided  M  ,  S, ,  and  S„  are  fixed. 
00   1       2 

For  small  perturbation  turrets,  a  similar  argument  can  be  made. 

V.   VISCOUS  FLOW  PHENOMENA  LEADING  TO  OPTICAL  DISTORTION 

A.   Types  of  Viscous  Flow  Phenomena 

Viscosity  dominates  and  determines,  at  least  in  part,  the  gross 
characteristics  of  several  different  types  of  flow.   These  include  boundary 
layers,  shear  layers,  open  port  cavities,  wakes,  and  shock  waves.   Although  not 
a  flow  of  the  same  nature  as  those  listed  above,  jet  engine  exhausts  should  be 
included;  exhausts  can  have  a  detrimental  effect  on  beam  propagation.   At  the 
speeds  and  geometric  scales  involved  with  aircraft,  the  viscous  flows  lead  to 
turbulence.   Extensive  literature  [6-21]  exists  on  the  propagation  of  electro- 
magnetic waves  through  a  turbulent  medium. 

Boundary  layers  surround  the  aircraft  surfaces.   The  boundary  layer  on 
leading  edges  and  forward  surfaces  is  laminar  with  transition  occurring  at  a 
critical  Reynolds  number.   Flow  separation  leads  to  shear  layers.   Shear  layers 
separate  two  flow  regions  with  different  flow  velocities. 

Due  to  the  extremely  high  irradiance  in  the  laser  beam,  the  turret  may  have  an 
open  port.   Windows  are  not  available  which  can  survive  the  high  irradiance.   The 
interior  of  the  turret  is  a  cavity  in  which  unsteady,  turbulent  flow  exists.   The 
flow  in  a  cavity  may  have  regularity  in  the  sense  that  discrete  vortices  can  be 
identified. 

Wakes  are  turbulent  and  adversely  affect  propagation.   When  selecting  a 
location  for  a  laser  turret,  any  location  within  or  near  a  wake  should  be  avoided. 


Shock  waves  create  large  changes  in  density  and  refractive  index.   Shock 
waves  which  originate  or  reflect  from  boundary  layers  create  intense  turbulence. 
Depending  on  conditions,  a  triple  point  may  occur  in  reflection  of  an  oblique 
shock.   Downstream  of  the  triple  point,  a  shear  layer  exists.   Also  downstream  of 
a  reflected  shock,  a  separated  flow  region  may  occur. 

Two  cases  can  be  identified.   In  one  case,  the  laser  beam  propagates  through 
a  region  disturbed  by  a  shock  wave.   In  the  other  case,  the  shock  wave  originates 
or  reflects  near  the  aperture  of  the  turret.   The  latter  case  causes  the  greater 
difficulty  with  optical  distortion. 

B.   Turbulent  Boundary  Layers 

The  boundary  layer  may  originate  on  aircraft  surfaces  far  upstream  of  the 
aperture.   If  this  is  the  case,  then  the  boundary  layer  will  be  thick.   When 
the  laser  turret  protrudes  into  the  flow,  the  boundary  layer  may  begin  on  the 
turret  itself.   In  that  case,  the  boundary  layer  is  thin.   One  can  state 
unequivocally  that  the  thinner  the  boundary  layer,  the  lesser  the  adverse 
influence  on  beam  quality. 

In  addition  to  the  distance  which  the  boundary  layer  has  traversed,  other 
factors  influence  the  boundary  layer.   Mach  number  characterizes  the  compressibility 
effects.   At  higher  Mach  numbers,  the  boundary  layer  tends  to  grow  in  thickness. 
Wall  curvature,  which  usually  is  accompanied  by  pressure  gradients,  has  an 
influence  on  boundary  layer  thickness  and  velocity  profile.   When  the  wall  is 
concave  to  the  external  flow,  Gortler  vortices  may  occur.   Cooling  or  heating  can 
alter  boundary  layer  properties.   In  fact,  cooling  may  be  one  approach  to  flow 
control.   Surface  roughness  may  amplify  turbulent  intensities.   Surface  roughness 
may  be  due  to  rivet  heads,  structural  joints,  insects  which  have  smeared  the  sur- 
face, etc. 
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Pressure  gradients  have  been  mentioned  earlier.   A  gradient  in  which  the 
pressure  increases  in  the  flow  direction  is  termed  an  adverse  gradient. 
Conversely,  a  gradient  with  decreasing  pressure  in  the  flow  direction  is  a 
favorable  gradient.   Adverse  pressure  gradients  lead  to  flow  separation  as  well 
as  thicker  boundary  layers.   Favorable  pressure  gradients  delay  transition  from 
laminar  to  turbulent  flow.   Favorable  pressure  gradients  tend  to  cause  thinner 
boundary  layers.   A  boundary  layer  in  a  favorable  pressure  gradient  rarely  separates 

Shock  wave,  boundary  layer  interaction  has  been  mentioned.   If  the  inter- 
action occurs  at  the  aperture  or  immediately  upstream  of  the  aperture,  severe  beam 
degradation  is  likely. 

For  turbulent  incompressible  boundary  layers,  the  thickness,  6,  is  given 
by  [22] 

-  =  0.37  Re  ~0,2  (22) 

x  x 

According  to  equation  (22),  6  increases  as  x  to  +  0.8  power.   For  a  laminar, 

compressible  boundary  layer,  the  thickness  is  given  by  [23] 

6  =  ClRex°'5  [1  +  C2  :LYJL  ^]  (23) 

where  C  and  C  are  constants.   Turbulent  compressible  boundary  layers  show 

similar  trends  with  increasing  Mach  number. 

A  boundary  layer  in  compressible  flow  has  a  temperature  gradient  across 

the  boundary  layer.   The  gradient  in  temperature  is  due  to  viscous  dissipation. 

A  quantity  known  as  the  recovery  factor  r  is  defined  for  an  adiabatic  wall;  r  is 

T  -  T 
r  =  TW  _    =  Pr  (24) 

0    °° 
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The  temperature  of  the  wall  is  T  ,  the  temperature  of  the  air  at  the  outer  edge 
of  the  boundary  layer  is  T^,  and  the  stagnation  temperature  is  T  .   Equation  (24) 
shows  r  equal  to  the  Prandtl  number,  Pr.   An  effective  Prandtl  number  can  be 
defined  for  a  turbulent  boundary  layer  and  will  have  a  value  near  unity. 
Hence  the  wall  temperature  will  be  equal  to  the  stagnation  temperature  when 
Pr  =  1.0. 

The  rms  density  fluctuations  in  the  turbulent  boundary  layer  are  related  to 

the  difference  in  density  at  the  wall,  p  ,  and  at  the  edge  of  the  boundary 

w 

layer,  p  .   According  to  Sutton  [24] 

p    =  0.1(p  -  p  )  (25) 

rms        °°    w 

for  boundary  layers  with  supersonic  external  flow.   For  subsonic  boundary  layers, 

the  constant  is  larger  than  0.1  shown  in  equation  (25).   The  density  at  the  wall 

is  given  by 

Poo 

pw  =  W~  (26) 

W 

The  rms  fluctuation  in  refractive  index  is  related  to  density  by 

n        p 
_rms  =  K^ms 

n        p  v 

where  k  is  the  constant  from  equation  (12). 

An  attenuation  coefficient,  a,  can  be  defined  for  propagation  through  a 
layer  of  a  turbulent  medium.   The  intensity  ratio  is 

f  -  e-°6  (28) 

0 

where  the  layer  thickness  is  6,  and  I   is  the  intensity  of  the  beam  entering  the 

layer.   The  attenuation  is  due  to  scattering  by  the  turbulent  boundary  layer. 

Sutton  [25]  gives  a  formula  for  a  as  follows: 
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2   2 
a  =  2k  n   A  (29) 

rms  v   ' 

where  k  is  the  wavenumber  equal  to  2tt/A  and  A  is  the  macroscale  of  the  turbulence, 
The  macroscale  is  obtained  from 


=  J    C(r)dr  (30) 


A 

0 

where  C(r)  is  the  correlation  function  for  refractive  index  in  an  isotropic 
turbulent  medium.   The  correlation  function  is  given  by 


/  An(r  +  R)An(r)47TR2dR 


C(r)  =  -2— (31) 


J     [An(R)] 


2   2 

4ttR  dR 


0 
where  An  is  the  deviation  from  the  average,  n, 

An  =  n(r)  -  n  (32) 

Sutton  [25]  reports  values  of  the  scale  of  turbulence  equal  to  about  0.1  6,  i.e., 
the  scale  of  turbulence  is  approximately  1/10  of  the  boundary  layer  thickness. 
For  propagation  through  turbulent  boundary  layers,  several  characteristic 
lengths  occur  as  follows: 

A     macroscale  of  turbulence 
6     boundary  layer  thickness 
D     aperture  for  laser  beam 
A     wavelength  of  laser  radiation 
Many  regions  for  turbulent  propagation  can  be  characterized  by  the  length  scales 
along  with  the  extinction  number  a6. 
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C.   Scaling  Propagation  Through  a  Turbulent  Boundary  Layer  of  Fixed  Thickness 

Two  parameters  can  be  used  to  characterize  the  degradation  in  beam  quality 

2 
due  to  turbulence.   One  is  the  decrease  of  intensity,  watts/cm  ,  on  the  axis 

of  the  beam  in  the  far-field.   For  large  values  of  D/A,  equation  (28)  provides  a 

means  to  calculate  I/In-   The  second  parameter  is  the  half-intensity  angle  for 

the  beam.   The  half-intensity  angle,  6-, /9>  defines  a  cone  with  a  vertex  angle 

equal  to  9  ,  .   On  the  surface  of  the  cone,  the  intensity  is  one-half  of  the  peak 

intensity. 

Figure  1  shows  the  half -intensity  angle  as  a  function  of  the  ratio  D/A.   The 
ordinate  is  the  variable  6  ,  /(1.129A/D) .   This  variable  compares  the  half- 
intensity  angle  for  the  scattered  beam  with  the  half-intensity  angle  for  a 
diffraction  limited  beam. 

Consider  a  C0?  laser  beam  at  10.6  microns  from  a  D  =  60  cm  mirror.   Assume  the 

boundary  layer  thickness,  6,  is  10  cm,  and  the  macroscale,  A,  is  1  cm.   For 

M  =  1.5  flow  with  turbulent  Prandtl  number  of  unity  and  p  =  1  atmosphere, 

p  =1.23  kg/m  ,  and  T  =  290°K,  the  values  of  p  and  p    are  0.848  kg/m  and 
00         &         °°  w      rms  & 

0.038  kg/m  .   Using  K  =  2.4  x  10~  ,  the  value  of  n   /n  is  7.44  x  10   .   The 

rms 

extinction  coefficient,  which  is  calculated  using  equation  (29),  is  0.39  per  meter. 
The  extinction  number,  a6,  is  0.039,  and  D/A  is  60.   The  corresponding  point  can 
be  plotted  in  Figure  1,  giving  6  ,  /(1.129A/D)  almost  on  the  dotted  line. 

Consider  now  an  iodine  laser  at  1.315  microns  with  a  turret  for  which 
D/A  =  S  has  been  maintained  constant.   Hence  D  =  7.44  cm;  also,  D/A  is  7.44. 
The  extinction  coefficient  grows  to 

a  =  (0.39)(8.06)2  =  25.34 
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TURBULENCE    CORRELATION  fuKCJjCN 

2"'       I     1,     >"»  •      3,    i 


SEPARATION    DISTANCE 
TURBULENCE  INTEGRA 


RATIO  OF  APERTURE    TO  INTEGRAL    SCALE  SIZE.  D/A 

Turbulence  corrcluliiui  function 


Figure  1.   Diffraction  angle  for  which  the  intensity 
is  one  half  the  maximum  value.   Reproduced 
from  Sutton  [25] . 


and  the  extinction  number  is  2.53.   The  point  for  the  iodine  laser  can  be  plotted 
in  Figure  1,  giving  a  half  angle  parameter  of  about  0.6.   The  increase  in  half 
angle  is  about  20  per  cent. 

For  a  thicker  turbulent  boundary  layer  with  6  =  20  cm,  the  results  are  more 
dramatic  as  shown  in  Table  II.   The  half  angle  parameter  increases  from  0.48  to 
1.2,  which  is  a  factor  of  2.5. 

It  should  be  emphasized  that  equation  (28)  does  not  apply  when  D/A  is  less 
than  approximately  10.   The  attenuation  is  less  than  that  predicted  by  equation  (28) 
when  D/A  <  10. 

D.   Scaling  Propagation  Through  a  Scaled  Turbulent  Boundary  Layer 

This  subsection  discusses  the  case  for  which  the  turbulent  boundary  layer 
originates  on  the  laser  turret  instead  of  an  aircraft  surface.   From  equation  (22), 
the  ratio  of  boundary  layer  thicknesses  is 
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Table  II.   Change  of  Half  Angle  with  Change  in  Wavelength 


wavelength,  A,  micron 

boundary  layer  thickness,  6,  cm 

macroscale,  A,  cm 

extinction  coefficient,  a,  per  meter 

extinction  number,  a.6 

ratio  D/A 

ratio  D/A 

Mach  number  M 

pressure,  p  ,  atmosphere 

3 
density,  p^,  kg/m 

o 
temperature,  T  ,   K 


10.6 

1.31. 

20 

20 

2 

2 

0.78 

50.67 

0.156 

10.13 

30 

3.72 

56,603 

56 

,603 

1.5 

1.5 

1 

1 

1.23 

1.23 

290 

290 

x  0.8  _   0.8 

6     x  L 

6  =  (/}    =  (L  } 

°2    Xl  Ll 


(33) 


Further,  the  ratio  of  macroscales  of  turbulence  is 


1  h. 

A1  =  61 


(34) 


The  extinction  number  scales  as 


n     ,2  .    T2   .   0.8  .   1.2 

a6)     A  A    L   L  L 

=  — -  — -  =  —  (—}  =    (—) 

a6)l    A22  Al   L22  V  (L2} 


(35) 


Equation  (35)  assumes  n    is  not  changed.   The  ratio  D/A  scales  as 
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^1   hh       h     h°'8  L2°'2 

D/A)1  "  4  A2  "  4  (L2}     =  (LX}  (36) 

Using  equations  (35)  and  (36)  along  with  data  of  Table  II,  the  following  values 
are  found  for  A  =  1.315  microns: 

D/A  =  19.8 

aS     =  1.89 
Plotting  the  preceding  points  in  Figure  1  yields  a  6  .   about  8  per  cent  larger. 
The  impact  of  reduced  A  is  much  less  when  the  boundary  layer  originates  on  the 
turret. 

VI.   SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

The  analysis  has  assumed  that  S  =  A/D  and  S_  =  L/D  were  maintained  constant 
as  A  was  decreased.   The  conclusions  are  valid  for  this  case. 

Using  the  optical  system  equation  of  Holmes  and  Avizonis  [3],  various  conse- 
quences of  reducing  A  were  examined.   In  summary,  the  following  should  be  noted: 

6       Phase  distortion  must  be  reduced  by  a  factor  of  8  for  a  change 
rms 

from  C0~  to  iodine  laser. 
K       A  factor  which  depends  on  the  spatial  distribution  of  intensity 

at  the  output  aperture  will  not  vary. 
T,      Diffractive  spillage  will  tend  to  increase  reducing  the  value  of 

V 

0       rms  jitter  can  be  maintained  constant  without  change  in  Iff- 
Detailed  comments  were  not  made  about  atmospheric  extinction  coefficient; 

however,  scattering  tends  to  increase  as  A  decreases. 

Turret  weight  and  volume  as  well  as  turret  aerodynamic  drag  will  be  reduced 

dramatically.   Much  greater  flexibility  in  the  number  and  location  of  turrets  is 

a  consequence.   Further,  the  possibility  of  obtaining  a  suitable  window  should 

improve;  open  port  problems  are  thereby  circumvented. 
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The  analysis  shows  that  the  phase  distortion  6/A  for  inviscid  flow  is  giver 
by 

y  =  (constant)  -r  (37) 

Since  D/A  is  constant,  the  phase  distortion  due  to  inviscid  flow  is  independent 
of  wavelength.   The  statement  is  true  provided  S  and  S   are  fixed.   Nonetheless 
severe  degradation  due  to  the  inviscid  flow  field  can  occur;  hence,  the  problem 
cannot  be  ignored. 

For  turbulent  boundary  layers,  the  value  of  the  extinction  number,  a6, 
increases  as  A  decreases.   Also,  the  ratio  D/A  decreases  as  A  decreases.   Since 
extinction  and  scattering  depend  on  a<5  and  D/A,  the  trends  are  not  favorable. 
The  half-intensity  angle  increases,  and  extinction  increases  as  A  decreases. 
Turbulent  boundary  layers  require  more  attention. 

Any  change  which  causes  a  thinner  turbulent  boundary  layer  tends  to  improve 
beam  quality.   Mounting  the  laser  turret  at  points  farther  forward  decreases  6. 
Cooling  the  wall  in  the  supersonic  case  also  should  be  an  improvement.   Finally, 
attention  to  surface  roughness  should  help  reduce  scattering  and  attenuation. 

When  the  turbulent  boundary  layer  originates  on  the  laser  turret,  the 
boundary  layer  thickness  decreases  as  the  turret  decreases  in  size.   Hence,  the 
impact  of  the  turbulent  boundary  layer  on  6. ,   is  almost  unchanged  by  a  change 
in  A;  however,  serious  degradation  in  intensity  due  to  a&   remains. 
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APPENDIX  A 
Consider  the  change  in  intensity  at  some  point  P  in  the  far-field  of  an 
aperture  illuminated  in  phase  and  amplitude  by 

F(£,n)  =  A(£,n)ejlK?'n)  (A-l) 

Silver  [26]  on  page  173  shows  the  amplitude  and  phase  of  the  wave  at 
point  P  is  given  by 


U  =  -1—  e 
P   XR  e 


"JkR  y  F(^n)exp[jk  sin  e(£  sin  <|)  +  n  cos  <$>)  ]d£dn     (A-2) 


A 
Intensity  at  point  P  is 


Ip  =  UpU*  (A-3) 


Introduce  nondimensional  variables 

£  =  S'D      n  =  n*D  (A-4) 

where  D  is  the  diameter  of  the  aperture.   The  integration  in  equation  (A-2) 

2 
extends  over  A  =  ttD  Ik.      Examine  first  the  exponent  in  the  integrand  of 

equation  (A-2);  it  becomes 

jk  sin  6[£  sin  (J)  +  n  cos  <$>]    =  jkD  sin  0[£'  sin  <J)  +  n '  cos  <p]  (A-5) 

Note  that 

kD  =  2tt/S1  (A-6) 

The  complex  function  F(£,ri)  is  exactly  the  same  as  equation  (A-l)  except  primes 

are  added  to  E,   and  n.   The  integral  in  equation  (A-2),  which  is  represented  by  I  , 

becomes 


D2In  =  D2  f   F(C,,n*)exp[j(27T/S1)sin  0(C  sin  $  +  n '  cos  <J>)]dC'dn'     (A-7) 
A' 

The  intensity  at  point  P  is  obtained  by  combining  equations  (A-2),  (A-3),  and  (A-7) 
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2  2     2 
D  I     I 

I   =  *  =  —5-  (A-8) 

P    \2   2  2  2 

Consider  two  cases  1  and  2  which  differ  in  the  wavelength  of  radiation  but  for 
which  S   is  constant.   The  ratio  of  intensities  becomes 

^  -  &  4  (a-9) 

P1    X„l  R2 
In  equation  (A-9) ,  the  integrals,  I  ,  cancel  when  P  is  located  at  the  same 
point  specified  by  R,6,<J).   Since  R  =  R  ,  the  ratio  of  intensities  is  identical. 
Consequently,  the  factor  K  in  equation  (1)  is  identical  for  cases  1  and  2. 
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